Role of isospin physics in supernova matter and neutron stars 
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We investigate the liquid-gas phase transition of hot protoneutron stars shortly after their birth 
following supernova explosion and the composition and structure of hyperon-rich (proto)neutron 
stars within a relativistic mean-field model where the nuclear symmetry energy has been constrained 
from the measured neutron skin thickness of finite nuclei. Light clusters are abundantly formed with 
increasing temperature well inside the neutrino-sphere for an uniform supernova matter. Liquid- 
gas phase transition is found to suppress the cluster yield within the coexistence phase as well as 
decrease considerably the neutron-proton asymmetry over a wide density range. We find symmetry 
energy has a modest effect on the boundaries and the critical temperature for the liquid-gas phase 
transition, and the composition depends more sensitively on the number of trapped neutrinos and 
temperature of the protoneutron star. The influence of hyperons in the dense interior of stars makes 
the overall equation of state soft. However, neutrino trapping distinctly delays the appearance 
of hyperons due to abundance of electrons. We also find that a softer symmetry energy further 
makes the onset of hyperon less favorable. The resulting structures of the (proto) neutron stars with 
hyperons and with liquid-gas phase transition are discussed. 



PACS numbers: 26.60.-c,21.65.-f 
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I. INTRODUCTION 

Recently considerable effort has been devoted to un- 
derstand the isospin physics because of its paramount 
importance in determining the structure and dynamics of 
finite nuclei [U-Q and also in the astrophysical contexts 
such as supernova dynamics, protoneutron star evolu- 
tions and the structure and properties of neutron star 
0-[3- I n stable nuclei, the neutron-proton asymmetry, 
S = (p n — Pp)/(p n + pp), is about 0.2, where p n and p p 
are the neutron and proton number densities. In forth- 
coming rare-isotope accelerator experiments the magni- 
tude of 5 could be well above 0.24. In contrast, a large 
S s=s 0.9 is reached in the dense interior of neutron stars 
USUI- Whereas, very small asymmetry of S ~ 0.1 is at- 
tained near the crust of (proto)neutron stars at densities 
p ~ lO^-lO" 1 fin" 3 In order to understand the 

composition and structure of stars over this extreme val- 
ues of density and asymmetry requires precise knowledge 
of the density dependence of nuclear symmetry energy 

^sym (p) • 

Several theoretical studies have demonstrated @, |1] 
that E sym (p) represents the leading coefficient of an ex- 
pansion of the total energy, E(p,S), for a neutron-rich 
system with respect to asymmetry: E(p, 6) = E(p, 0) + 
E sy m{p)S 2 + 0(5 4 ). The symmetric nuclear matter equa- 
tion of state, E(p, 0), has been already well constrained 
over a wide density range from analysis of a large body of 
data within non-relativistic models [Til [l2j and relativis- 
tic models [l[. In contrast, though the nuclear symmetry 
energy at the normal nuclear matter density po = 0.15 
fm~ 3 is well constrained to E sym (po) = 32 ± 4 MeV, its 
value and trend at lower and especially at higher densities 
are poorly known (H. [l3|. V arious theoretical calculations 
based on microscopic 14j4l6|| an d/or ph enomenological 
many-body approaches [2, U, [HI, [T34l3 > though repro- 



duce the well-constrained value of E sym (po), their high 
density predictions of symmetry energy are extremely 
diverse. While certain models predict a monotonically 
increasing symmetry energy with density, other model 
results show that E sym (p) increases initially up to about 
Po and decreases thereafter. These diverse predictions 
in the density dependence in the symmetry energy lead 
to large uncertainty in the incompressibility of symmet- 
ric nuclear matter which mainly stems due to insufficient 
knowledge about the isospin dependence of in-medium 
nuclear effective interactions. 

Recently some progress has been achieved by 
consistently constraining the symmetry energy to 
32(p/p ) ' 7 < E sym (p) < 32(p/p ) 1 ' 1 at subsaturation 
densities from analysis of isospin diffusion [H, [2(| and 
isoscaling [2lj data in intermediate energy heavy ion colli- 
sions and from the study of neutron skin of several nuclei 
22. 12311 . Within a relativistic mean-field (RMF) model 
17.l24j[25|. th e sy mmetry energy has been constrained to 
a small range [23j from analysis of neutron skin thickness 
of several nuclei across the periodic table. Further, the 
slope parameter and the incompressibility for the con- 
strained E sym (p) was found [23] consistent with that ex- 
tracted from isospin diffusion and isoscaling data. 

In an effort to pindown the precise density dependence 
of symmetry energy, it has been noted [3, [26| that the 
skin thickness of neutron-rich heavy nucleus, for example 
208 Pb, are linearly correlated with the symmetry pres- 
sure at a density slightly below the normal nuclear matter 
value. Since the same pressure supports a neutron star at 
supranormal densities against gravitational collapse, sug- 
gests an underlying correlation between neutron skin of 
heavy nuclei and neutron star radii [26[ . Within the RMF 
models, the constrained E sym (p) demonstrated [H[ that 
the neutron skin in 208 Pb tend to yield smaller neutron 
star radii. Furthermore, within the constrained range of 
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symmetry energy, significant sensitivit y o n several fea- 
tures of liquid-gas phase transition |27l. |28| ] in hot asym- 
metric nuclear matter was found 29]. In particular, the 
boundary and area of the liquid-gas coexistence region, 
the maximal isospin asymmetry and the critical values of 
pressure and isospin asymmetry all of which systemati- 
cally increase with increasing overall softness of E sym (p). 

It is thus instructive to explore the effects of symme- 
try energy in the RMF model, that has been limited to 
a narrow range from analysis of skin data, on the com- 
position and structure of hot protoneutron and cold neu- 
tron star matter that encompasses a large density range. 
A newly born neutron star or protoneutron star (PNS) 
[H, IH, [13, H3| is formed in the aftermath of a successful su- 
pernova explosion and subsequent gravitational collapse 
of the core of a massive star. Due to short mean free 
path of neutrinos of A„ ~ 10 cm (compared to PNS radii 
of about 10 km), the neutrinos are trapped temporarily 
within the neutrino-sphere of the star. The PNS are thus 
rich in leptons, mostly e~ and v ei and quite hot with cen- 
tral temperatures of T = 10—50 MeV. Thus the star mat- 
ter has the interesting possibility to undergo liquid-gas 
phase transition close to the surface of the star [3l|, [32| , 
the characteristics of which is mainly controlled by the 
symmetry energy at subsaturation densities. In the much 
longer Kelvin-Hclmholtz evolution stage of about 10 sec, 
the neutrinos escape from the interior, and the hot and 
lepton-rich PNS changes into a cold and deleptonized 
neutron star. The loss of neutrinos enforce the electrons 
and protons to combine resulting in a neutron-rich neu- 
tron star matter. Consequently, the symmetry energy at 
supranormal densities plays an important role in deter- 
mining the composition, in particular the possible pres- 
ence of strange hyperons, the structure (mass and radius) 
and the long-term cooling of the neutron star. 

In this paper we study the role of isospin on the liquid- 
gas phase transition in hot lepton-rich protoneutron star 
and the properties of (proto)neutron star within the ac- 



curately calibrated parameter set FSUGold 25[ in the 
relativistic mean-field (RMF) model. For this purpose, 
we restrict the model density dependence of symmetry 
energy E sym (p) within a narrow range that reproduced 
the skin thickness of several nuclei [231 ]. At the supra- 
normal densities the contribution of the hyperon degrees 
of freedom to the equation of state of the star is con- 
sidered. At very low densities well below the saturation 
value, that is relevant to the crust, few-body correlations 
become important. The system can then minimize its 
ener gy v ia the formation of medium-modified light clus- 
ters [33L H3] which eventually dissolve into homogeneous 
nuclear matter at higher densities due to Pauli blocking. 

The paper is organized as follows: In Sec. II we intro- 
duce the RMF model for baryons along with the medium- 
modified clusters. The construction of the liquid-gas 
phase transition in the hot protoneutron star matter is 
presented. In Sec. Ill the numerical results for the struc- 
ture and composition of the star matter are discussed. 
Sec. IV summarizes the results. 

II. FORMALISM 



The RMF Model with clusters 



We consider the relativistic mean-field (RMF) model 
with all the charge states of the baryon octet (B = 
p, n, A, E + , S°, S~, S°, S _ ) and their antiparticles. Light 
clusters (C = 2 H, 3 H, 3 He, 4 He = d,t,h,a), that 
are treated as quasi-particles, are also included as ex- 
plicit degrees of freedom. The properties of these clus- 
ters are modified by medium effects via the temperature- 
dependent shifts in the binding energies as treated in Ref. 
[33L l3~ij . The interaction Lagrangian density in the non- 
linear RMF model [25|, [U, Hfl for a system of baryons, 
clusters and leptons for protoneutron star matter is given 
by 



i=B,t,h 

-\™f€*^ + \(d^a mW) - I (g aN af - A { g aN(J f - + l - m ^^ + | (gl^^f 

-\p»v pf 1 " + \m 2 pP ^ p» + A v (g 2 pN p^ p»){gl N u^) + ( V ~ ™*) (1) 

1 

I 



where ipi is spin-1/2 fields for the particles (i = B,t,h), 
cj) a is the spin-0 field for the alpha particle, and </>^ is the 
spin-1 field for deuteron all of which are coupled to the 
<J,cu,p mesons. The field strength tensors for the vector 
mesons are u)^ v = d^ L uj v - <9„w M and p M „ = <9 M p„ - d v p^. 
The nonlinear cr-meson couplings (k, A) soften the sym- 



metric nuclear matter equation of state (EOS) at around 
Po, while its high density part is softened by the self- 
interactions (£) for the w-meson field. The sum on I is 
over the non-interacting electrons and muons (e~ and 
p~ ) and their antiparticles in the star. 
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The covariant derivative 



iff 1 - guiAiu" - gpilNi - Z t \r- f^hi , (2) 



of the particle i with mass, proton, neutron numbers 
(Ai, Z%, Ni) involves interaction with the vector mesons. 
Of course for hyperons (i = Y) the (Ai, Zi, Ni) drop out 
in Eq. @ and hereafter. I 3i is the isospin projection of 
the baryon and cluster charge state i. The effective mass 
of the ith species is 



g ai AiG - AB,, 



(3) 



where the vacuum rest mass of the clusters i = d,t,h,a is 
given by m, = Zim p + A^m„ — Bf. For the vacuum bind- 
ing energy Bf of the clusters we adopt the experimentally 
measured values. Note the clusters which are predomi- 
nant at densities p ~ (10 — 10 )po are treated as struc- 
tureless point particles. The scalar isovector ao(980) (the 
8 meson) that splits the effective masses of the baryons 
has not been included in the present study as its effects 
on the star properties were found to be minimal [H [3(| . 
The medium dependent shift of the binding energy ABi 
that appears only for the clusters are adopted from Ref. 
34]. For nuclei embedded in nuclear matter, the shift has 
been calculated by solving the in-medium Schrodinger 
equation which leads to the ground state energy [34J 



Ba,z — B Al 



K 2 
2Ar, 



AB 



SE 
A,Z 



AB 



Pau 
A y Z ' 



AB 



Cou 
A,Z- 



(4) 



The shift due to self-energy, AB^ Z , is already contained 
in the effective particle mass in the RMF model. The 
Coulomb shift, AB C A °%, that can be evaluated within the 
Wigner-Seitz approximation gives a small value for the 
light clusters and thus neglected. The dominant con- 
tribution emerge from the Pauli shift, AB F A a % , that has 
been evaluated in the perturbation theory with Jastrow 
and Gaussian approaches for d and (t, h, a), respectively. 
We adopt the empirical form for Ref. [3J] 



tot n tot 
r 



-Pi 



1 



mm 



SBt(T), (5) 



where p z = 2{Z l p t ° t + Np^/A,, and p°(T) = 
B$/6Bi(T). The quantity 5Bi(T) can be found from 
Eqs. (26)- (27) in Ref. [Hj]. With increasing temperature 
and/or density a cluster will dissolve when its total bind- 
ing energy Bi — B® + ABi vanishes. In the above expres- 
sion the total densities for neutrons, p^ ot = J2i=B c ^»P*> 
and protons, Pp 0t = Yli=B C %iPi are contributions from 
baryons (B) and clusters (C). 

The composition of the star is determined by the re- 
quirements of charge neutrality and /3-equilibrium con- 
ditions B\ — > Bi + I + 77; , B 2 + I — > Bi + vi , where Bi 
and Bi are the baryons and I is a lepton (e - or pr) 
or their respective antiparticles at finite temperature. 
With clusters (A, Z), these weak processes correspond to 
(A,Z) -> (A,Z+l)+l+Vt and (A,Z)+l -> [A,Z-l)+vi. 



The chemical potential of a species i with baryon num- 
ber Bi, charge Qi, and lepton number L t in chemical 
equilibrium is in general 



Mi 



BiHB + QiPQ + Lip L . 



(6) 



The three independent chemical potentials pb, pq and 
Pl can be determined from total baryon number, charge, 
and lepton number conservations leading to 



PA.z = AflB + ZflQ, 

Pe- = -p e + = -p Q + p L , 

p Vl = -p-t = PL] (l = e~,fT). 



(7) 



For a deleptonized neutron star, the (anti)neutrinos es- 
cape freely so that pl = 0. The remaining chemical po- 
tentials, ps and pq, can be obtained from conservations 
of total baryon and charge densities 



tot i ~tot 
Pv ' 



p = A ^ = Pn 
i=B,C 

PQ= + 

i—B^C l—e,fi 



Qwi = 0. 



(8) 



When neutrinos are trapped, the matter (inside the 
neutrino-sphere) on a dynamical time scale is character- 
ized by a fixed lepton fraction (number of leptons per 
baryon), Yli = Yl + Y Vl , for each flavor I = e, p. In 
fact, muons are absent in a PNS that leads to Yl^ = 
and p Vl = p Ve in Eq. ([7| . Consequently, the composi- 
tion of a hot PNS is characterized by a constant value of 
Ylc — Y e + Y Ue , which in the late stage reaches a value 
of -0.4. 

The number density for a fermionic species i — B,t,h,l 
at temperature T^0 



Pi 



2J, + 1 
(2tt) 3 



d 3 k [f?(k)-fr(k)], 



(9) 



where Ji is the spin of the ith species and the distribution 
function for particles and antiparticles (referred to as ± 
sign) is as usual [2J], |28| 



ft(k) = [exp{(E*(k)Tv i )/T}+l}- 1 



(10) 



The density of the bosonic species i — d, a is obtained 
from 34 1 



2J. + 1 
(2tt) 3 



d 3 k gi{k) 



(11) 



where pi stems from the condensed bosons in the ground 
state. The Bose-Einstein distribution follows as 



g i (k) = [exp{(E*(k)-v i )/T}-iy 



(12) 



For interacting baryons and clusters, the effective en- 
ergy and chemical potentials are respectively E* = 
^Jk 2 + m* 2 and v { = p A ,z -g^A.^- g pl (Ni- Zi)p 03 I 3i , 



4 



where ujq and po3 are the timelike and isospin three com- 
ponent of the vector w-meson and vector-isovector p- 
meson. For field-free lepton, these quantities reduce to 
E\ = Ei = ^k 2 + m 2 and V\ — pi. The field equations 
for the mesons (m = a,&, p) are derived in the usual way 
for the Lagrangian of Eq. (JTJ) . Herein the contribu- 
tion from the cluster densities to the source terms in the 
meson fields are consistently incorporated by noting that 
in the present Lagrangian the couplings g m i of the me- 
son fields to the baryons and clusters (i = B ,d,t,h,a) 
are constant (density independent) unlike that in Ref. 
0. 

At finite temperature and baryon density, the energy 
density and pressure corresponding to the Lagrangian of 
Eq. ([T]) can be obtained from the thermodynamical po- 
tential D, 27] as 



£ 



P 



\ 22 
— — + y {g*Na) + ^ (g*N<j) + 

+ E 2 4^- [d 3 kEt(k)[f+(k) + fr(k)] 



-B,t,h,i 



E 

i—d,a 



(2tt) s 

2 J, + 1 
(2tt) 3 



d 3 k E*(k) gi {k) +m*pi 



(13) 



m 2 a a 2 k 3 A 4 m 2 w 2 
2 ^| V9<?N<J) - (gaNcr) + 

+ 7^ (g^N^o) 4 H + A„ (gojN^af (g P NPm) 2 

+ E l^/ d3fc i^[/^) + /r(fc)] 



=B.t,h. 



3(2tt) 



E f§^ / 



A: 2 



3i(A;). 



(14) 



where x = V G /V being the fraction of the total vol- 
ume occupied by the gas phase. The stability condition, 
Eq. (fTS)) , implies the inequalities [27J d 2 T /dp 2 > and 
(d 2 F/dp 2 )(8 2 T/dp 2 3 ) > (d 2 F/dp 3 dp) 2 . These inequali- 
ties are equivalent to 



dP 



(d 2 T 



dp ) T s 9 \dp 2 



dpp 
85 



T,S 

< or 



> 0, 



T,S 



T,P 



dp n 

88 



> 0. 



(17) 
(18) 



T,P 



The first inequality indicates mechanical stability which 
means a system at positive isothermal compressibility re- 
mains stable at all densities. The second inequality cor- 
responds to chemical stability which suggests that energy 
is required to alter the concentration in a stable system 
while maintaining temperature and pressure fixed. When 
one of these conditions get violated, a two-phase system is 
energetically favorable. The coexistence liquid-gas phase 
is governed by the Gibbs criteria for equal pressures and 
chemical potentials in the two phases with different den- 
sities but at the same temperature: 



P(T,p L ) = P(T,p G ) 
^{T,p L ) = p t (T,p G ) 



(19) 
(20) 



where the chemical equilibrium condition refers to all the 
species (i — n,p,d,t,h,a,e,v e ) in the system. Since we 
are dealing with free leptons, it is evident from Eq. (|2"0)l 
that the lepton density in the liquid and gas phases are 
identical and thus forms a common background in the 
baryonic matter. Else for interacting leptons, pf ^ p G , 
and thus the lepton number conservation would result in 
an additional condition similar to those of Eq. (TTF 



C. Symmetry energy and model parameters 



B. Liquid-gas phase transition in PNS matter 

The RMF model is applied to study the liquid-phase 
(LGP) transition in hot protoneutron star matter. Such 
a coexistence phase may occur at subsaturation densities 
where the star is composed of nucleons ( N ) , clusters (C) 
and leptons only. A single-phase system at temperature 
T and density p is stable against LGP separation if its 
free energy density J 7 is lower than the coexistence liquid 
(L) and gas (G) phase: 



T(T,p,p 3 ) < (1 



X )T L {T,p L ,p L 3 ) 



The total baryon density, p = p P ° 
density, p 3 = p\ 
satisfy 



tot 
in 



+ X J 7G (T,p G ,p^). 

(15) 

9^ ot , and the isospin 
-p5, for a N-C-e-v e system 



P = (1 — x)p + XP 

P3 = (1- X)P3 +Xpf, 



(16) 



The original FSUGold [25[ , with an isovector coupling 
Ay = 0.03, produces a soft symmetric nuclear EOS with 
an incompressibility Kq = 230 MeV. In addition to the 
ground state properties of several nuclei, this model re- 
produces the strengths of the giant monopole resonance 
in 90 Zr and 208 Pb. The nucleonic matter symmetry en- 



ergy for this Lagrangian at temperature T 
by 



Esym(p) 



k 2 



9pN k 



FN 



6E 



FN 



12tt 2 



is given 



(21) 



FN 



FN 



m^ 2 , and the Fermi momen- 



where E 

turn and effective mass for nucleons are fci?jv an d Tn* N = 
m N — go- n<^, respectively. The effective p- meson mass 
is m* 2 = m 2 + 2g 2 N (A v g 2 N ujQ). In asymmetric nuclear 
matter, the density dependence of symmetry energy has 
been studied [23L HH by simultaneously varying A„ and 
g P N so that for all combinations of (A„, g p w) the sym- 
metry energy is fixed at E sym (j5,T = 0) = 26 MeV at 
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an average density ~p corresponding to kppf = 1.15 fm _1 . 
Note the symmetry energy even at the saturation density 
has not been accurately constrained experimentally and 
only an average E syia at po is constrained by the binding 
ener gy o f nuclei. Thus by following the above prescrip- 
tion [26[ one ensures accurate binding energies and pro- 
ton densities of heavy nuclei such as 208 Pb. In fact, only 
those A v values were selected that reproduced the two 
measured observables mentioned above. For the present 
study we adopt A v = 0.0 — 0.03 since the corresponding 
E sym (p,T = 0) and their slopes and curvatures are in 
reasonable agreement with that extracted from neutron 
skin thickness of several nuclei as well as the isoscaling 
and isospin diffusion data 23[. It may be noted that 
with increasing A„ the density dependence of symmetry 
energy becomes overall softer, i.e., softer at supranor- 
mal densities and stiffer at p < po- Within this limited 
A„ range, considerable effect on the features of liquid-gas 
phase transition was also observed in hot asymmetric nu- 
clear matter (29j . 

The hyperon- meson coupling constants, g a y and g U Y, 
are strongly correlated by the potential depth of the hy- 
peron Y in saturated nuclear matter Q as 



U< 



(AT) 



(22) 



We use SU(6) symmetry for the vector coupling constants 



1 



1 



1 




p(fm ) 



P(fm ) 



:3wS — 5wE 



FIG. 1: (Color online) Density dependence of particle frac- 
tions Yi at temperatures T = 1,2,4,6,10,20 MeV in sym- 
metric nuclear matter (SNM) and in protoneutron star (PNS) 
matter at a fixed lepton fraction of Ylb = 0.1 in the FSUGold 
model with coupling = 0.03. 



9 P N = ^ffps = Sph; 9 P A = 



(23) 



The g a Y couplings are adjusted to reproduce the esti- 
mated hypernuclear potential depth Uy at po. From 
detailed studies of AN interaction, an attractive poten- 
tial depth of ?7 A = — 30 MeV has been estimated [13] so 
that the observed bound A hypernuclear states are repro- 
duced. Several S hypernuclear states [H| and free Ss pro- 
duced in reactions |39l . |40| suggest an attractive poten- 
tial of ui N) = -18 MeV. Though the A and S potential 
depths have been well studied, however, data on E hyper- 
nuclei are scarce and ambiguous due to strong EiV — >• AN 
decay. While previous studies of hyperon stars have used 



u 



(AT) _ 



-30 



an attractive E potential (with Uf~ 
MeV) that are based on the observed bound |He hy- 
pernucleus, recent analysis 0, HU of E~ atomic data 
point to a sizable repulsive nuclear E potential of depth 



u. 



(AT) 



-30 MeV. For the present study we adopt this 



repulsive depth of the E potential. This ambiguity will 
modify considerably the hypernuclear composition in the 
star, which however, has smaller effects on the bulk prop- 
erties of /3-equilibratcd star matter. 



III. RESULTS AND DISCUSSIONS 

We begin this section by presenting in Fig. [T] the 
fractions, Y, = AiPi/p, of free protons and clusters 



i = p,d, i, h, a at various temperatures for symmetric 
nuclear matter (SNM: left panels) and neutrino-trapped 
protoneutron star (PNS: right panels) matter in the 
FSUGold relativistic mean-field model for the isovector 
coupling A v = 0.03. The overall trend exhibited by sym- 
metric matter is quite similar to that for neutron-rich 
protoneutron star matter. At very low densities p < 10~ 7 
where the meson fields and the baryon chemical poten- 
tials drop rapidly the PNS behaves as SNM with Y p = 0.5 
(see also Fig. |3J). With the formation of the clusters, 
the threshold density of which is quite sensitive to the 
system temperature, the free proton fraction decreases 
rapidly. At a given temperature, the two-nucleon cor- 
relation results in the formation of deuteron first, fol- 
lowed by three-particle h, t and finally the a particle. 
At T < 7 MeV, the strongly bound alpha is seen to be 
the most abundant particle over a wide density range 
as many-body correlations gain prominence with density. 
Moreover, with increasing temperature, the a fraction de- 
creases while other cluster abundances increase; the peak 
positions gradually shifting to higher densities. Conse- 
quently with increasing T, the minimum of Y p increases 
and shifts to higher densities and then decreases again 
at T > 7 MeV when deuterons dominate. This was also 
observed in the RMF model with density-dependent cou- 
plings |34j. Note that at a given temperature and den- 
sity the (free) neutron-rich supernova matter has some- 
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what smaller cluster fractions compared to the symmet- 
ric nuclear matter. The clusters eventually dissolve due 
to Pauli blocking at higher densities p > 10~ 2 fm~ 3 re- 
sulting in a homogeneous matter with Y p = 0.5 for SNM 
(left panel) and a neutron-rich n-p-e-v PNS matter (right 
panel) . 



16 - 



12 - 



4 - 



= 0.00 



1 1 

^ (w/o clusters) 



\ 
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FIG. 2: (Color online) Boundary of liquid-gas phase transition 
in a supernova matter at a fixed lepton fraction of Vie = 0.4 
in the FSUGold model with couplings A„ = 0.0 and 0.03 with 
and without the inclusion of clusters. The inset shows the 
symmetry energy effects for A v = 0.0 and 0.03 at the liquid- 
phase boundary without clusters. 

We now consider the possibility of liquid-gas phase 
(LGP) transition in hot and neutrino-trapped protoneu- 
tron star consisting of neutrons, protons, electrons and 
electron neutrino. Figure [2] shows the boundary for the 
existence of the coexistence phase at various tempera- 
tures in the FSUGold model for the isovector coupling 
A„ = 0.0,0.03. The results are for electron lepton frac- 
tion Ylc = Yl = 0.4 (and muon lepton number Yj, M = 0) 
which approximately corresponds to the value at the on- 
set of trapping. At a nucleon density of p w 0.08 fm -3 , 
the critical temperature for A„ = 0.03 is T c = 14.6 MeV 
which is close to that for symmetric nuclear matter [29| . 
Whereas a softer nuclear symmetry energy E sym at sub- 
saturation density (A„ = 0), gives a somewhat smaller 
T c . In fact, the critical temperature divide the coex- 
istence boundary into a low-density gas phase (on the 
left) and a high-density liquid phase (on the right) and 
the liquid-gas phase coexist within the boundary. Due 
to trapping, the net electron concentration is quite high, 
and the matter at the low baryon density phase is con- 
trolled by non-degenerate electrons. As a consequence, 
the nuclear symmetry energy corresponding to A„ = 
and 0.03 has practically no effect on the gas phase bound- 
ary. On the other hand, due to dominant contribution of 
the degenerate nucleons in the relatively high-density liq- 
uid phase, larger effects from symmetric pressure in the 
stiffer E Bym for A„ = 0.03, shifts the liquid-phase bound- 
ary to a slightly higher density (see inset of Fig. [5]). The 
protoneutron star matter is found to exhibit a LGP tran- 
sition even at very low densities at a high temperature of 
T ~ 1 MeV. 



With the inclusion of light clusters (solid line in Fig. 
[5]), the softening of the EOS shifts the low-density coexis- 
tence phase boundary to a higher value at temperatures 
T > 4 MeV where clusters become abundant. On the 
other hand, the high-density boundary that is dominated 
by free nucleons is relatively unaffected. Thus the clus- 
ters effectively reduce the width of the phase coexistence 
region and increases the critical temperature to T c = 15.7 
MeV for A„ = 0.03 case. 
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FIG. 3: (Color online) Composition of supernova matter with 
trapped neutrinos at a lepton fraction of Yl £ = 0.4 at a 
temperature of T = 1 and 2 MeV in the FSUGold set for 
Ai, = 0.03 coupling. The results are for uniform matter (solid 
line) and with liquid-gas phase transition (dashed line). 

In Fig. |3J the particle fractions are shown for neutrino 
trapped matter at an electron lepton fraction Y^ e = 0.4 
at temperatures of T = 1 and 2 MeV for uniform mat- 
ter (solid line) and in the coexistence liquid-gas phase 
(dashed line). One of the major effects of trapping is 
to enhance the net electron fraction Y e over the entire 
baryon density range; the magnitude of Y e gradually in- 
creases with Yl- Considering first uniform matter (solid 
line), at p < 10 -7 fm -3 , the energy and pressure of the 
hot protoneutron star matter is entirely dominated by 
the non-degenerate electrons because of their small mass. 
Due to charge neutrality, the rapid rise of net electrons 
at this low density enforces an equally large enhancement 
of net number of free protons and a consequent drop in 
free neutron density due to baryon number conservation 
(see Figs. [3ja) and[3][c)). This tends to make the lepton- 
rich matter in this density regime highly symmetric with 
identical free neutron and proton fractions. 

For uniform matter at higher densities, clusters dom- 
inate the particle yield (compared to free protons and 
neutrons especially at high T; (Figs. [3Jb) and[3Jd)) so 
that the charge neutrality is maintained primarily by the 
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bound protons in these clusters. Compared to Y^ e = 0.1 
(Fig. [TJ), larger electron fraction in the present = 0.4 
case enhances the cluster production via charge neutral- 
ity. At all temperatures of interest in the supernova dy- 
namics, we find a and deuteron to be the dominant clus- 
ters. Thus, within this density range the star matter is 
also highly symmetric but due to abundant symmetric 
cluster formation. 
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FIG. 4: (Color online) Density dependence of nuclear asym- 
metry S for supernova matter with trapped neutrinos at var- 
ious lepton fraction of Yt e at a temperature of T = 1 MeV. 
The results are in the FSUGold set at A„ = (0.0, 0.03) for 
uniform matter (dashed, solid line) and with liquid-gas phase 
transition (dotted, dash-dotted line). 

With the onset of liquid-gas phase transition, it is im- 
portant to note from Fig. [3] that cluster formation drops 
drastically (dashed lines). This stems from vanishing 
cluster contribution from the liquid part of the mixed 
phase that has a density well above the transition value 
where the clusters dissolve, as well as due to small clus- 
ter fraction in the gas phase that has a density smaller 
compared to the total p. Thus at small T < 4 MeV 
in the coexistence phase, Y p is effectively enhanced and 
becomes identical to Y e . This leads to the formation of 
n-p-e-v matter that tends to be rather symmetric even 
in the coexistence phase. However, with increasing T, 
as the cluster fractions are large and the onset of phase 
transition shifts to higher p (see Fig. [5]), considerable 
amount of clusters can occur in the low-density uniform 
phase of a hot protoneutron star. The formation of clus- 
ters and near symmetrization of the PNS matter has two 
major consequences. First, during the evolution of the 
supernova matter, the abundant light nuclei produced 
within the neutrino-sphere 0, [To[ at densities p < 10~ 5 
fm~ 3 may be ejected outside and contribute to the nu- 
cleosynthesis of heavy elements. Second, the decrease in 
the pressure due to the nuclear symmetry energy over- 
whelms the total increase from the thermal and leptonic 
pressures. This makes the equation of state softer with or 



without phase transition compared to the neutrino-free 
matter. 

It may be noted from Fig. [3l that at a fixed Yl, anti- 
neutrino D e dominates over its particle counterpart at 
lower densities and/or in the coexistence region. More- 
over, the trapping density is also shifted to a higher value 
in the coexistence phase. Thus, apart from neutrino 
transport, the exact location of the neutrino-sphere is 
quite sensitive to the influence of the coexistence phase 
in the supernova matter. 




io" 7 io" 5 io" 3 io" 1 io" 6 io" 4 10" 2 10° 
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FIG. 5: (Color online) Per particle energy E/A, free energy 
F/A, and pressure P/p as a function of density for supernova 
matter at a lepton fraction Yhe = 0.1 and at temperatures 
T = 0.5,1,2 MeV. The results are in the FSUGold set at 
couplings A v — 0.0 (solid lines) and 0.03 (dotted lines) for 
uniform matter and with liquid-gas phase transition. 

Figure. [4l shows the variation of neutron-proton asym- 
metry 5 = (/o^ ot — Pp 0t )/p with baryon density at several 
fixed lepton fraction Yi, e that correspond to the delep- 
tonization epoch of the supernova at T — 1 MeV. Con- 
sidering the A„ = 0.03 set, the dip seen in the asymme- 
try for uniform matter (solid lines) at p ~ 10~ 2 feu 
is due to cluster formation that increases with Y^ e . A 
further rapid drop in the asymmetry at p % 10 fm -3 
results from lepton dominance in the star. Thus, even in 
the late stage of neutrino diffusion (at small Yl), abun- 
dant nuclei can be produced in the near-symmetric mat- 
ter [3l], HH . With the onset of liquid-gas phase transition 
at pth ~ 0.1 fm -3 (dash-dotted line), the asymmetry is 
found to decrease faster compared to the uniform matter. 
This stems from dominant contribution from the liquid 
part of the mixed phase which exhibits a dropping asym- 
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metry at p > pth (solid lines) plus a fraction from the 
very low-density region of the gas phase that is nearly- 
symmetric. Both these effects systematically favor a free 
energy with small 8. Thus, compared to uniform mat- 
ter, the supernova matter at T ~ 1 MeV is found to be 
nearly symmetric in the coexistence phase at larger den- 
sities near the neutrino-sphere and the PNS behaves as 
n-p-e-v matter. Of course at higher temperatures, the 
cluster contribution would be enhanced and the onset of 
coexistence phase is also delayed to a higher density; see 
Figs. [3fc)-(d). At a given p, this causes the PNS to 
be somewhat less symmetric compared to that at lower 
temperatures. 

The sensitivity of E sym on asymmetry 6 has been also 
studied by considering the A v = set which gives an 
overall stiffer density dependence of symmetry energy. 
In this case, for uniform matter (dashed line), a larger 
(smaller) 6 is obtained at densities smaller (larger) than 
the normal nuclear matter value. However, the effects 
of symmetry energy is negligible both at small densities 
that are dominated by leptons, and also in the mixed 
phase (dotted line). 

Figure [5] depicts the equation of state (EOS) of su- 
pernova matter at Y^ e = 0.1 for uniform matter and in 
the coexistence phase for A^ = 0.0 (solid lines) and 0.03 
(dotted lines) with clusters. At supranormal densities, 
the matter is characterized by repulsive w-meson inter- 
action that causes a rapid increase in the EOS. At about 
the normal density of po ~ 0.15 fm -3 , the uniform mat- 
ter shows a minimum for relatively small temperatures. 
This is a manifestation of the saturation properties of nu- 
clear matter when thermal effects contribute modestly 
A larger E sym for the A„ = 0.03 at p < po (dotted line) 
gives a smaller binding energy. In fact, the sensitivity of 
symmetry energy on the EOS persists till the contribu- 
tion from free nucleons and/or clusters dominate up to a 
density of p ~ 10 -4 fm -3 . 

In the coexistence phase, the (free) energy drops ap- 
preciably with decreasing density at small temperatures. 
As mentioned above, this is a consequence of the liquid 
part that contributes to enhance the binding of the sys- 
tem. The overall qualitative features found in the ther- 
modynamic quantities with clusters are similar to those 
without clusters for the supernova matter (not shown). 
Compared to nucleons only supernova matter, the inclu- 
sion of cluster degrees of freedom reduces somewhat the 
total free energy. Further, the total pressure is found 
to be small at p ~ 10~ 5 fm~ 3 when most of the clus- 
ter species contribute collectively. The pressure eventu- 
ally becomes stiffer when only the massive alpha con- 
tributes dominantly at p ~ 10~ 2 fm~ 3 (see inset of Fig. 
[7]) . Consequently, clusters delay the onset of the coexis- 
tence phase (situated in the soft low p regime of the EOS) 
to a higher density compared to that without clusters. 

We now consider the constitution and structure of 
beta-equilibrated and charge neutral matter in the su- 
perdense regime. It may be noted that thermal effects 
contribute negligibly in this density regime where the nu- 




P (fm" 3 ) p (fm" 3 ) 



FIG. 6: (Color online) Relative concentration of beta- 
equilibrated matter containing nucleons and hyperons in cold 
lepton-poor matter (Y v = 0) and lepton-rich matter (Yl b = 
0.4) at finite entropy per baryon S/B = 1 in the FSUGold set 
at At, = 0.0 and 0.03 values for the coupling constant. 



cleons become degenerate. In Fig. [SJ the abundances of 
baryons and leptons as a function of density are shown 
for the FSUGold parameter sets A„ = and 0.03 with 
and without the trapped neutrinos. The composition 
of the matter should be quite sensitive to the symme- 
try energy E sym as it relates to the chemical potential 
via fi = fx n — fj, p ~ 4E sym S 0, [|[ which in turn, deter- 
mines the hyperon production threshold density by the 
condition ps = — 1b (Pe ~ Pv e ) > £fl, where eb is 
the energy density of the baryon species B. For neutron 
star matter with Y v = 0, the stiffer density dependence 
of symmetry energy at A„ = causes larger enhance- 
ment of electron (and hence proton) fraction compared 
to A v — 0.03 case. However, with the appearance of 
negatively charged S - hyperon which competes with lep- 
tons in maintaining charge neutrality, the lepton (e - and 
/i~) concentrations begin to fall. Consequently, the sen- 
sitivity of E sym on the composition is diminished. The 
charge neutral S° hyperon then appears, which occurs at 
a somewhat earlier density for the stiff symmetry energy 
compared to the soft A„ = 0.03 set. 

Trapped neutrinos have a large influence on the compo- 
sition of the protoneutron star as evident from Fig. [5] for 
Yle = 0.4. The appearance of hyperons is now delayed to 
higher densities since in the above mentioned threshold 
condition, the term (/i e — fi Ve ) turns out to be positive. 
With the emergence of A hyperon, the dropping of proton 
and thereby the electron fraction is now compensated by 
the rise of v e abundances so that Yl c is maintained con- 
stant. A softer symmetry energy (A„ = 0.03) is found to 
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FIG. 7: (Color online) Pressure versus energy density for f3- 
equilibrated matter with free nucleons and clusters only (npc) 
and with further inclusion of hyperons (npcH). The results are 
for cold lepton-poor matter (Y v = 0) and lepton-rich matter 
(Vie = 0.4) at finite entropy per baryon S/B = 1 in the 
FSUGold set at A„ = 0.0 and 0.03 values for the coupling 
constant. The inset shows the equation of state of the crust 
for uniform matter (dashed line) and with liquid-gas mixed 
phase (solid line) at A„ = 0.03. 



favor a significantly larger Y Ve . Incidentally, at very high 
densities (fx e — fx Ve ) becomes negative. The impact of this 
alteration in the lepton chemical potentials is the emer- 
gence of positively charged E + before the E°'~ hyperon. 
From the central densities of the maximum mass stars 
(see Table I), it is evident that for the repulsive potential 

potential of = +30 MeV used here, E hyperons will 
not be present in the star matter at all. On the other 
hand, if we adopt an attractive interaction potential of 



E/j: = —30 MeV, as used in previous studies, the E~ 
hyperon will appear at a density of (2 — 3)po followed by 
E° and E + in the star matter. In such a situation Ss will 
not appear in neutrino-trapped matter. 

The equation of state (EOS), i.e. pressure P versus en- 
ergy density e is displayed in Fig. [7] for the neutrino- free 
matter (Y v ) and matter with an electron lepton fraction 
Ylb = 0.4 corresponding to that of Fig. [51 For nucleons- 
only star, trapped neutrinos makes the EOS softer as 
the decrease in the symmetry pressure due to enhanced 
proton fraction exceeds the increase in leptonic pressure. 
Compared to the stiff A„ = set, a smaller symmetry 
energy /pressure at high density for the A„ = 0.03 set ob- 
viously produces a softening in the EOS. However, as a 
soft E sym (p) favors a large v e fraction [3(|, the softening 
is reduced in the lepton-rich matter. With the inclusion 
of hyperons, the EOS becomes softer in both lepton-rich 
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FIG. 8: (Color online) The mass-radius relation for star se- 
quences for matter with free nucleons and clusters only (npc: 
solid line) and with further inclusion of hyperons (npcH: 
dashed line) in a cold lepton-poor matter (Y v — 0) and lepton- 
rich matter (Yl = 0.4) at finite entropy per baryon S/B — 1 
in the FSUGold set at A v — 0.0 and 0.03 values for the cou- 
pling constant with BPS equation of state [45[ at the crust. 
At Yl = 0.4, the results for npc star at A„ = 0.03 with 
an uniform matter (dotted line) and liquid-gas mixed phase 
(dash-dotted line) EOS at the surface are also shown; see text 
for details. 



and lepton-poor matter. This is essentially caused by the 
release of the Fermi pressure of the nucleons into the hy- 
peronic degrees of freedom, and also due to decrease of 
pressure exerted by the leptons as these are replaced by 
negatively charged hyperons [5] . However in a protoneu- 
tron star, the trapped-neutrinos inhibit the onset of hy- 
perons (see Fig. [5]) which makes the EOS with hyperons 
stiffer compared to the neutrino-free matter — a rever- 
sal of behavior relative to the EOS for nucleons-only star 
matter. As evident, the comparative stiffness increases 
for the larger E synl (A v = 0) that allows relatively slightly 
larger hyperon population in the cold catalyzed neutron 
star. 

The differences in the EOS should be reflected in the 
structure of the stars, namely their masses and radii. 
While most of the mass originate from the dense inte- 
rior of the star, the crust has a negligible contribution 
of ~ 10 _5 Af Q to the total mass. In contrast, about 
40% of the radius originates from the EOS at p < po 
and thus should be strongly influenced by the possible 
liquid-gas phase transition at this density regime. At 
first, the star sequence is calculated with the low density 
equation of state taken from Baym-Pethick-Sutherland 
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TABLE I: Maximum masses, M max /M0, their corresponding radii, RM max , an d their central densities, p c /po, of stars without 
(Y v = 0) and with (Yz,e = 0.2, 0.4) trapped neutrinos; the normal nuclear matter density is po = 0.148 fm . Results are in the 
FSUGold mean-field model at A„ = and 0.03. 
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Y Le 


= 0.2 


7.082 


1.749 


11.322 


6.067 


1.444 


12.083 




Y Le 


= 0.4 


7.420 


1.737 


11.108 


6.716 


1.569 


11.674 




Yv ~- 


= 


7.815 


1.722 


10.840 


5.814 


1.344 


11.953 


0.03 


Y Le 


= 0.2 


7.815 


1.711 


10.747 


6.856 


1.459 


11.341 




Y L e 


= 0.4 


7.618 


1.710 


10.952 


6.858 


1.578 


11.502 



(BPS) of Ref. [isf : the effects of the coexistence and pure 
phase with clusters are discussed later. Here, the over- 
all EOS is constructed by smoothly matching (without 
any discontinuity in pressure, baryon and energy den- 
sities) the BPS equation of state at p < 0.001 fm -3 , 
followed by the mixed phase EOS comprising of leptons 
and free nucleons only for the inner crust of the star in 
the RMF model, and finally with the RMF model EOS 
at p > po, with or without hyperons. The static star se- 
quence obtained by solving Tolman-Oppenheimer-Volkoff 
equations 43j are shown in Fig. [5] for the different equa- 
tions of state. Listed in Table I are the maximum masses 
M max , corresponding radii i?M max , and the central bary- 
onic densities p c . With free nucleons and leptons-only 
star, the trapped-neutrinos give smaller maximum mass 
M max (compared to Y v = case) as the pressure from 
the stiff EOS can support larger masses against gravita- 
tional collapse. With the introduction of hyperons, the 
reversal of behavior observed in the EOS, i.e. trapping 
leads to a stiffer EOS, is manifested in larger maximum 
mass protoneutron stars. Note in general, the inclusion 
of hyperons leads to a relatively smaller mass stars due to 
softening of the EOS. The larger radii in stars with hyper- 
ons are a consequence of weaker gravitational attraction 
from smaller masses that causes the stars to be large and 
diffuse. For the neutrino-free stars, though A„ = 0.03 
gives smaller masses than A„ = 0, the corresponding star 
radii are however found to be smaller in the former case. 

The baryonic interactions in this model raises the max- 
imum mass from M max w O.7M for the EOS of a free 
gas [43| by a factor of about two for nucleons-only stars. 
However, the inclusion of hyperons in cold neutrino-free 
stars yield maximum masses that are smaller than the 
current observational lower limit of IAAMq imposed by 
the larger mass of the binary pulsar PSR 1913 + 16 [44] . 

The inset of Fig. [7] shows the EOS for hot matter 
with entropy per baryon S/B = 1 and Yl c =QA for uni- 
form matter (dashed line) and in the coexistence liquid- 
gas mixed phase (solid line). For clarity, the EOSs cor- 
responding only to the high-density (liquid) regime of 
the mixed phase at volume fraction \ = V G /V > 10~ 4 
(see Eq. (|15p ) are shown. In fact, \ can be as low as 



10~ 8 near the low-density (gas) phase boundary. Figure 
[8] illustrates the effect of liquid-gas phase transition at 
low density on the mass-radius sequence of free nucle- 
ons and clusters-only protoneutron stars for A.„ = 0.03. 
Though the maximum masses are identical with M max w 
1.710Mq the radii in the sequence of stars are all dif- 
ferent. In particular, a stiff EOS obtained in the mixed 
phase from Gibbs construction gives consistently large 
radii protoneutron stars. Compared to the BPS [45[ 
equation of state used above at p < 0.001 fm -3 , the 
RMF models in general provide much stiffer EOSs at the 
crust that result in stars with much larger radii. 



IV. SUMMARY AND CONCLUSION 

In this article, within an accurately calibrated relativis- 
tic mean field model where the nuclear symmetry energy 
has been constrained from measurements of neutron skin 
thickness of finite nuclei, the structure and properties of 
neutron and protoneutron stars are studied. For hot and 
lepton-rich protoneutron stars we examine the possible 
liquid-gas phase transition near the crust comprising of 
free nucleons and light clusters. Compared to free nu- 
cleon abundances, light clusters are found to dominate 
the particle yield at moderate and high temperatures in 
an uniform supernova matter. Further, the fraction of 
trapped neutrino is shown to play an important role in 
making the supernova matter more symmetric via en- 
hanced production of deuteron and alpha clusters. With 
liquid-gas phase transition, the yields of clusters drop 
drastically at temperatures T < 4 MeV in the coexis- 
tence phase; here near symmetrization over a wide den- 
sity p ~ 10~ 5 — 10 _1 fm~ 3 is achieved by the free nucleons 
that are not bound in the clusters. At these low densities, 
symmetry energy has modest impact on the coexistence 
phase boundaries and on the equations of state of mat- 
ter compared to thermal effects and on the number of 
trapped neutrinos. 

The influence of symmetry energy on hyperon produc- 
tion in the dense interior of cold neutron stars and hot 
lepton-rich protoneutron stars are studied. Within the 
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present knowledge of the potential depth of hyperons in 
bulk nuclear matter, we find a repulsive potential of X 
hyperons inhibits its appearance in the star matter. Neu- 
trino trapping delays the appearance of hyperons due 
to relative disposition of chemical potentials of various 
species. A softer symmetry energy in general result in 
smaller chemical potentials that delays further the onset 
of hyperons and thereby reduces the hyperonization of 
the star matter. The emergence of hyperons soften the 
nuclear equation of state considerably in the FSUGold 



model. This results in stars with very small masses es- 
pecially for the cold deleptonized matter. We also find, 
that compared to uniform matter, a stiffer equation of 
state in the coexistence liquid-gas phase near the crust of 
hot lepton-rich matter yield stars with larger radii. The 
ongoing and future hypernuclear programs can provide 
accurate information of the potential depth of hyperons 
in bulk nuclear matter, which in conjunction with bet- 
ter mass and radius measurements of neutron stars could 
constrain the nuclear symmetry energy more accurately. 
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